Available  online  at  www.sciencedirect.com 


ELSEVIER 


%V  ScienceDirect 


Journal  of  Power  Sources  161  (2006)  949-956 


JOURNAL  OF 


www.elsevier.com /locate /jpowsour 


Effects  of  cathode  humidification  on  the  gas-liquid 
interface  location  in  a  PEM  fuel  cell 

Chun-I  Lee*  1 ,  Hsin-Sen  Chu  * 

Department  of  Mechanical  Engineering,  National  Chiao  Tung  University,  Hsinchu  300,  Taiwan,  ROC 

Received  25  March  2006;  received  in  revised  form  5  May  2006;  accepted  8  May  2006 

Available  online  1 1  July  2006 


Abstract 

This  study  investigates  the  location  of  the  gas-liquid  interface  under  various  humidity  conditions  in  the  cathode  gas  diffusion  layer  and  the 
conventional  flow  field  of  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC),  using  commercial  computational  fluid  dynamics  (CFDRC)  analysis, 
based  on  a  three-dimensional  steady  state  isothermal  model.  The  effects  of  the  formation  of  liquid  water  on  the  transport  of  reactant  gas  are 
considered  in  the  model.  The  simulation  results  indicate  that  the  gas-liquid  interface  approaches  the  gas  flow  channel  inlet  and  the  cell  performance 
declines  gradually  as  liquid  water  obstructs  the  pores  in  the  porous  media  as  the  relative  humidity  of  the  cathode  is  increased.  Additionally,  the 
oxygen  fraction,  the  water  fraction  and  the  liquid  water  saturation  field  along  the  flow  channel  are  presented. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  promis¬ 
ing  new  energy  conversion  devices  because  of  their  high  effi¬ 
ciency,  simplicity  of  design  and  operation  and  environmental 
friendliness.  However,  the  performance  and  efficiency  of  this 
PEM  fuel  cell  must  be  improved  further  to  make  it  cost-effective 
for  use  in  these  applications  [1-3].  Accordingly,  water  man¬ 
agement  is  critical  to  the  proper  operation  of  PEM  fuel  cells. 
However,  over  the  past  decade,  effective  water  management 
has  remained  elusive.  During  operation,  formed  liquid  water 
is  transported  within  the  membrane-electrode  assembly,  and 
controlling  its  motion  is  extremely  difficult.  When  the  rate  of 
generation  of  water  at  the  cathode  by  electro-osmotic  drag  and 
the  oxygen  reduction  reaction  exceeds  the  water  removal  rate 
from  the  cathode  by  back-diffusion  to  the  anode,  evaporation, 
and  capillary  transport  of  liquid  water  through  the  cathode  dif¬ 
fusion  layer  and  catalyst  layer,  the  cathode  becomes  flooded. 
Additionally,  excess  accumulated  liquid  water  blocks  the  gas 
pores  in  the  diffusion  layer  and  the  catalyst  layer,  which  are 
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needed  for  the  transport  of  oxygen  gas,  and  form  a  barrier  over 
the  catalyst  active  surface  in  the  catalyst  layer,  worsening  the 
performance. 

Springer  et  al.  [4,5]  developed  a  ID  isothermal  model  to  pre¬ 
dict  the  net  water  per  proton  flux  ratio  across  the  membrane  in 
PEMFCs.  However,  this  model  cannot  evaluate  the  decrease  in 
reactants  and  the  accumulation  of  products  in  the  flow  direction. 
Fuller  et  al.  [6]  and  Nguyen  et  al.  [7]  proposed  a  2D  heat  trans¬ 
port  and  water  transport  model  to  account  for  the  distribution 
of  water  and  the  net  water  flux  within  the  membrane,  at  vari¬ 
ous  operating  temperatures  and  membrane  hydration  along  the 
flow  channels.  Yi  et  al.  [8]  modified  Nguyen’s  model  to  describe 
both  liquid  and  gas  phase  distributions  along  the  flow  paths  on 
the  anode-side  and  the  cathode-side  in  a  PEMFC.  The  results 
indicated  that  anode  humidification  improved  the  conductivity 
of  the  membrane,  and  that  the  cell  performance  was  enhanced  by 
the  injection  of  liquid  water  and  the  increase  in  humidification 
temperature. 

In  Refs.  [9,10],  Mazumder  et  al.  developed  a  3D  model  to 
predict  the  effect  of  the  formation  of  liquid  water  on  the  per¬ 
formance  of  the  PEMFC.  In  Ref.  [9],  the  polarization  behavior 
was  overestimated,  based  on  the  assumption  that  no  liquid  water 
was  formed.  In  contrast,  when  liquid  water  was  present,  the 
condensation  and/or  evaporation  of  water  via  equilibrium  phase 
transformation  proceeded  as  described  in  Ref.  [10].  Since  water 
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Nomenclature 

C 

mass  fraction 

CF 

quadratic  drag  factor 

D 

mass  diffusivity  (m2  s-1) 

F 

Faraday  constant  (C  mol-1) 

I 

current  density  (A  m-2) 

kc 

condensation  rate  constant  (s-1) 

kt 

evaporation  rate  constant  (atm-1  s-1) 

kp 

permeability  (m2) 

M 

molecular  weight  (kg  mol-1) 

P 

pressure  (atm) 

R 

universal  gas  constant  (J  mol-1  K-1) 

S 

source  term 

T 

temperature  (K) 

u 

velocity  in  the  v-direction  (ms-1) 

V 

velocity  in  the  y-direction  (ms-1) 

V 

operating  voltage  (V) 

w 

velocity  in  the  z-direction  (ms-1) 

xw 

molar  fraction  of  water 

Zf 

spices  valence 

Greek  letters 

a 

transfer  coefficient  for  the  reaction 

£ 

porosity 

X 

water  content  of  membrane 

P 

density  (kg  cm-3) 

a 

ionic  conductivity  of  the  ionomer  (£2-1  m-1) 

X 

tortuosity  of  diffusion  layer 

V 

viscosity  of  flow  (kg  m-1  s-1) 

0 

phase  potential  (V) 

Subscripts 

a 

anode 

c 

cathode 

channel 

flow  channel 

CL 

catalyst  layer 

eff 

effective 

GDL 

gas  diffusion  layer 

in 

inlet 

k 

kth  component  of  fuel  reactant 

Mem;  m 

membrane 

sat 

saturated 

X-direction 

y 

T-direction 

z 

Z-direction 

Superscripts 

a 

anode 

c 

cathode 

ref 

reference 

condensation  can  obstruct  pores  and  reduce  cell  performance, 
the  polarization  behavior  of  PEMFC  must  be  accurately  pre¬ 
dicted  to  prevent  such  a  situation.  Restated,  the  consideration  of 
liquid  water  transport  substantially  enhances  the  ability  to  pre¬ 


dict  cell  performance,  and  improves  the  consistency  between 
the  experimental  data  and  simulation  results,  especially  at  high 
current  densities. 

Wang  et  al.  [1 1]  proposed  a  multiphase  mixing  model  to  study 
two-phase  water  transport  in  the  cathode.  Their  results  revealed 
that  the  formation  of  liquid  water  influences  the  performance  of 
the  cell.  Nevertheless,  the  overall  PEMFC  simulations  could  not 
be  understood  further  because  they  considered  only  the  cathode- 
side.  You  et  al.  [12]  presented  a  2D  two-phase  flow  mixture 
model  to  analyze  the  flow  and  transport  in  the  cathode  of  the 
PEMFC,  and  to  predict  the  phase  change  and  the  water/oxygen 
distribution  in  both  the  flow  channels  and  the  porous  diffusion 
layer.  Also,  such  a  two-phase  model  provides  actual  characteris¬ 
tics  of  the  flow  and  the  transport  of  reaction  gas  for  thermal/water 
management. 

More  recently,  computational  fluid  dynamics  (CFD)  have  uti¬ 
lized  water-transport  models  to  better  predict  cell  performance 
at  high  loads.  Yan  et  al.  [13,14]  presented  a  2D  numerical  model 
for  the  reactant  gas  transport  phenomena  and  cell  performance 
with  baffle  effects  in  the  flow  channel  of  the  bipolar  plate.  The 
effects  of  the  baffle  width  are  explored.  Additionally,  the  effects 
of  liquid  water  formation  on  the  reactant  gas  transport  are  taken 
into  account  in  the  modeling  and  examined  in  the  analysis. 

Recently,  most  of  the  research  on  the  humidification  con¬ 
ditions  has  focused  on  cell  performance  [8,15,16].  Studies  in 
this  area  reveal  that  the  relative  humidity  critically  affects  cell 
performance.  As  pointed  out  by  Natarajan  et  al.  [16],  the  per¬ 
formance  of  the  cell  increases  as  the  relative  humidity  of  the 
cathode  decreases. 

This  study  investigates  the  location  of  the  interface  when  the 
liquid  water  began  to  condense  along  the  flow  channel  of  PEM- 
FCs;  its  contribution  to  the  cell  performance  is  also  discussed. 
Additionally,  the  resulting  oxygen  and  water  fraction  distribu¬ 
tions  and  liquid  water  saturation  fields  at  fixed  cathode  humidity 
were  obtained  to  validate  the  simulation  results.  All  equations 
were  solved  using  the  commercial  computational  fluid  dynamics 
(CFDRC)  software  to  simulate  this  system. 

2.  3D  mathematical  model  of  PEMFC 

Fig.  1  schematically  depicts  the  3D  PEMFC  system  employed 
in  this  study.  It  consists  of  the  anode  flow  channel,  the  anode 
diffusion  layer,  the  anode  catalyst  layer,  the  proton  exchange 
membrane,  the  cathode  catalyst  layer,  the  cathode  diffusion  layer 
and  the  cathode  flow  channel.  Table  1  presents  the  geometri¬ 
cal  and  physical  parameters  adopted  in  the  present  simulation; 
some  minor  parameters  are  omitted  and  can  be  found  elsewhere 
[10,18,19].  In  this  study,  the  inlet  fuel,  hydrogen,  is  assumed  to 
be  fully  humidified  on  the  anode- side;  the  operating  pressure  is 
1  atm  and  the  operating  temperature  is  343  K.  The  objective  is 
to  investigate  the  location  of  the  interface  when  the  liquid  water 
begins  to  condense  along  the  flow  channel  of  PEMFCs.  The  fol¬ 
lowing  simplifying  assumptions  are  made:  (1)  the  gaseous  phase 
of  the  working  fluid  behaves  as  an  ideal  gas  and  the  liquid  water 
is  incompressible;  (2)  the  Reynolds  number  of  the  fluid  is  below 
100  because  the  velocity  of  the  mixture  is  low,  and  the  flow  is 
considered  to  be  laminar;  (3)  the  diffusion  layer,  the  catalyst 
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Fig.  1.  Schematic  diagram  of  3D  PEMFC. 


Table  1 

Geometric  and  physical  parameters  in  the  simulation 


Channel  length  (m) 

0.07112 

Channel  height  (m) 

7.62  x  10“4 

Channel  width  (m) 

7.62  x  10“4 

Diffuser  layer  thickness  (m) 

3  x  1(T4 

Catalyst  layer  thickness  (m) 

1  x  10“5 

Membrane  thickness  (m) 

3  x  10“5 

Porosity  of  diffuser  layer 

0.4 

Permeability  of  diffuser  layer  (m2) 

1.76  x  1CT11 

Porosity  of  catalyst  layer 

0.4 

Permeability  of  catalyst  layer  (m2) 

1.76  x  1CT11 

Porosity  of  membrane 

0.28 

Permeability  of  membrane  (m2) 

i.i8  x  icr18 

Operation  temperature  (K) 

343 

Operation  pressure  (atm) 

1 

Anode  fuel 

Hydrogen 

Cathode  fuel 

Air 

Hydrogen  stoichiometric  flow  rate 

1.5 

Air  stoichiometric  flow  rate 

3 

Relative  humidity  of  the  anode  (%) 

100 

layer  and  the  membrane  are  isotropic  and  homogeneous,  and 
characterized  by  an  effective  permeability  and  uniform  poros¬ 
ity;  (4)  the  system  operates  in  a  steady  and  isothermal  state;  (5) 
negligible  heat  is  generated. 

Based  on  these  assumptions,  the  basic  gas  transport  equations 
for  a  general  3D  PEMFC  are  as  follows: 


Continuity  equation: 


du  dv  dw 

- 1 - H - =0 

dx  dy  dz 

Momentum  equations: 


£eff 


du  du  du 
u  - - b  v - b  w 


dx  dy 
sQ  ff  dP 


dz 


p  dx 


d2u  d2u  d2u 
+  Ueeffl  dx2  +  dy2  +  dz2 


+  Su 


£eff 


dv  dv  dv 
u  - — b  v - b  w 


dx  dy 

£eff  dP 


dz 


P  fy 


d2v  d2v  d2v 
+  USeffl  dx2  +  dl2  +  d? 


+  5 


V 


£eff 


dw  dw  dw 
u - b  v - b  w 


dx  dy 
£eff  dP 


p  dz 


+  V£q  ff 


dz 
d2W 


d2W  d2W 

dx2  dy2  dz2 


Species  concentration  equation: 
dCk  dCh  dCk 

£eff  l  M— - h  V  — - b  w 


dx 


dy 


dz 


n  .  3 2Ck  d2Ck  d2Ck 
=  Okfi{{  |  + 


+  S 


w 


dx2  dy: 


dz2 


Sc  Sl 


Table  2  presents  expressions  for  the  source  terms  Su,  Sv,  Sw  and 
Sc  in  the  above  equations.  Within  the  table,  the  parameters  eeff» 
Cp,  kp  and  Zf  are  the  effective  porosity,  the  quadratic  drag  factor, 
the  permeability  and  the  valence  of  the  species,  respectively; 


Table  2 

Expressions  for  the  source  terms  in  the  governing  equations 


’w 


Flow  channel 
Diffuser  layer 

Catalyst  layer 


Membrane 


0 


/Cp 


/Cp 


V£ 


i^L^/u2  +  v2  +  w2 

Cl«y?l^u2  +  V2  +  W2 

sS 


- r —  U  ~ 


effvpu 


Vw2  +  V2  +  w2  + 


v 


Zf  Cjp-  F  •  V  (p 


du 

dx 


0 


-f-w 

/Cp 


-f-w 

/Cp 


/Cp 


fj=p\nr  +  w2  +  to2 

fk^V»2  +  u2  +  w2 


S^Vw2  +  l>2  +  U,2  + 


V^p 

2Zf  Ch+F-V0 


dv 

dy 


o 


-gw-  ^SV»2  +  v2  +  to2 


\fT’ 


-Fw~ 

/Cp 

e5sElfl^u2  +  v2  +  w2 

-gw- 

/Cp 

2  +  v2  +  w2  + 


P 


b  ZfCH+  F  V<P  ■  *ijjr 


0 


^2  :  2  FCa^a’  ^2  : 

H20  I  2FC~cJc 


ZF  n 

AT  ^C,eff,H+ 


CH+ 


92</>  ,  a2</>  i  d2^ 


a^2 


a>-2 


3z2 
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Dk  e ff  =  DkSTi  represents  the  effective  diffusion  coefficient  of 
the  Mi  component  of  fuel  reactant  [17] ;  ja  and  jc  are  defined  as: 


( 


ChA 


e(«a  F/RT)r] 


1 

e(ac  F/RT)r] 


( 


cbA 


Q(aaF/RT)r] 


1 

e(ac  F/RT)rj 


The  phase  potential  equation  yields  the  local  current  density 
distributions: 


d  f  d&\  d  f  d  f  d&\ 

to  Vm  9^ )  +  dy  Vm7 )  +  3z  )  =  Sj  (8) 


where  0  and  crm  represent  the  phase  potential  and  the  ionic 
conductivity  of  the  membrane. 

The  relationships  between  the  phase  potential  0  and  the  cur¬ 
rent  density  i  can  be  expressed  as, 


ix  =  -a. 


m 


d0 

dx 


ly  -  CTm 


d0 

dy 


lz  -  <Jm 


d0 

3  z 


(9) 


Then,  Eq.  (8)  can  be  rewritten  as 


3  ix  3  iy  3  iz 

dx  dy  d z  1 


(10) 


In  the  membrane,  the  ionic  conductivity  crm  is  strongly  related 
to  the  water  content  A.  It  is  defined  as  the  ratio  of  the  number  of 
water  molecules  to  the  number  of  charge  sites  [4] : 


<?m(T)  =  o+  exp 


(11) 


+ef  =  0.005 1  39a  -  0.00326 


(12) 


The  empirical  relationship  between  the  water  content  in  the 
membrane  and  the  partial  pressure  of  the  water  is, 


0.043  +  17.81a  —  39.85a2  +  36.0a3  0  <  a  <  1 

“  “  (13) 

14  +  1.4(a  —  1)  1  <  a  <  3 


where  a  is  the  water  activity  and  is  given  by, 


(14) 


In  the  above  equation,  the  saturation  pressure  varies  with  the 
temperature  and  can  be  determined  from  the  thermodynamic 
table  or  using  the  following  empirical  expression: 


jq— 2.1794+0.029537 — 9.1837x  10_572  +  1.4454x  10_773 


(15) 


During  the  operation  of  the  fuel  cell,  the  partial  pressure  of  the 
water  in  the  electrode  may  exceed  its  saturation  pressure  if  the 
local  concentration  of  water  is  high.  Accordingly,  liquid  water 
may  form  and  occupy  the  pores  in  the  porous  media.  Operating 
the  cell  at  a  high  reaction  rate  may  cause  severe  mass  transport 
overpotential  because  the  diffusing  species  are  blocked.  Further¬ 
more,  extremely  small  pores  in  the  porous  media  cause  capillary 
forces  to  dominate  the  transport  of  liquid  water.  However,  the 
actual  expression  of  this  force  cannot  be  formulated.  Eq.  (16) 


models  the  equilibrium  phase  change.  This  approach  has  also 
been  adopted  by  Nguyen  et  al.  [17],  and  is  the  same  approach 
being  used  in  current  study,  because  it  allows  natural  extension 
to  the  finite-rate  phase  change. 

The  determination  of  liquid  water  condensation  or  evapora¬ 
tion  capacity  can  be  denoted  as  Sfi  [10]: 


Mkc 


£eff-^w 

RT 


(*w  P 


kQ£Qffsp(xw  P  Psa  t), 


ifvwP>  Psat 
ifvwP  <  Psat 


(16) 


where  kc  and  ke  are  the  condensation  and  evaporation  rate  con¬ 
stants,  respectively,  s  is  the  saturation  of  liquid  water,  and  is 
defined  as  the  ratio  of  liquid  volume  to  the  volume  of  the  pores. 
vw  P  is  the  partial  pressure  of  water  vapor,  and  vw  is  the  molar 
fraction  of  water  vapor.  The  saturation  pressure  of  water  is  Psat, 
and  can  be  computed  using  the  curve-fitted  expressions  provided 
by  Springer  et  al.  [4]. 


2.7.  Boundary  conditions 

Wang  et  al.  [11]  indicated  that  although  the  gas  in  the  GDL 
is  transported  by  diffusion,  some  convection  with  a  very  low 
velocity  also  occurs.  Therefore,  the  boundary  conditions  in  this 
work  are  expressed  as  general  expressions.  In  this  work,  the 
operating  temperature  and  pressure  are  set  to  343  K  and  1  atm, 
respectively.  A  constant  flow  configuration  was  used  with  gas 
fluxes  of  1 .5  and  3  times  the  stoichiometric  value  at  1  A  cm-2 
for  H2  and  air,  respectively. 

At  the  gas  channel  inlet, 

anode  :  u  =  0,  v  =  0,  w  =  iufn,  Ch2  =  Ch2  in’ 

Ch2o  =  C&2  0,in  (17) 


cathode  :  u  =  0,  v  =  0,  w  =  wfn,  Co2  =  Cq2  in, 

6+t20  —  6+i20,in’  6"n2  =  6+j2,in  (^) 

At  the  interfaces  between  the  channel  walls  and  plate  collectors, 

3  Ck 


u  =  V  =  w  = 


dx 


=  0 


(19) 


At  the  interfaces  between  the  gas  channels  and  the  gas  diffuser 
layer, 


Af+GDL 


+ff,GDL 


+ff,GDL 


3  u 

3  u 

dy 

y=^GDL 

dy 

dv 

dv 

dy 

y=^GDL 

dy 

3  w 

3  w 

y — Ehannel 


y — Ehannel 


dy 


y=^GDL 


dy 


Uy=Y GDL  —  Uy= ^channel 


y — Ehannel 


7 — tGDL  ^ y — Ehannel  ’  OJy — Fqdl  0)y — ^channel  ’ 


-eff.GDL 


3 Ck 
dy 


y=^GDL 


3 Ck 
dy 


y — Ehanncl 


^k,y=YQDL  ~  Ck,y=Ychmne\ 


(20) 
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At  the  interfaces  between  the  gas  diffuser  layer  and  the  catalyst 
layer, 


The  boundary  conditions  of  phase  potential  at  the  interface 
between  the  catalyst  layer  and  the  membrane  is, 


£eff,GDL 


£eff,GDL 


£eff,GDL 


du 

dy 

dv 

dy 

dw 

dy 


du 

=  £eff,CL  — 

y=FGDL 

dv 

=  £eff,CL  — 

y=rGDL  dy 

dw 

=  Self, CL  — 

y=rG  DL 


}-=Ycl 


y=YCL 


WV=Egdl  uy=Yc  L’ 


y=Ycv 

vy= Egdl  =  vy=YcL’  wy=Yc,DL  =  wy=  Ecl  » 


£eff,GDL 


dCk 

dy 


dCk 

=  £eff,CL 

y=FGDL  dy 


y=rcL 


C/c,y=YaDL  Clc,y=Yc L 


(21) 


At  the  interfaces  between  the  catalyst  layer  and  the  membrane, 


£eff,CL 


£eff,CL 


du 

dy 

dw 

dy 


dv 

=  0,  ^eff,CL  — 

y=ECL  °y  y=YCi, 

0,  Uy=Y(2  L  ^5 


=  0, 


y=Ycv 


vy=YCL  =  Wy=YCL  =  0,  £eff,CL 


3Q 

dy 


=  0  (22) 


;y=tcL 


The  effect  of  liquid  water  upon  cell  performance  is  investigated. 
Proton  transport  and  the  electro-osmotic  phenomenon  cause 
liquid  water  to  be  transported  in  the  proton  exchange  mem¬ 
brane.  Therefore,  in  elucidating  the  transport  of  liquid  water, 
the  velocity,  the  flux  and  the  shear  stress  at  the  catalyst  layer 
and  the  membrane  interface  must  be  considered  to  be  continu¬ 
ously  distributed.  Additionally,  the  saturation  concentration  of 
liquid  water,  Cs,  and  the  concentration  flux  must  be  continuous. 
Therefore,  the  following  boundary  conditions  were  applied. 
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vy=YwiEM  —  vy=YcL  ’  Emem  —  wv=Ycl 


^eff,MEM 


dCs 

dy 


dC  s 

=  £eff,CL  — 
J=Emem  y 


y=ECL 


^S,y=EMEM  ^S,V=Ecl 


(23) 


At  the  gas  channel  outlet,  the  fully  developed  flow  conditions 
are  assumed  to  be 


^cl  =  ^MEM, 


£eff,  CL 


d@CL 

dy 


9^mem 

dy 


(25) 


3.  Results  and  discussion 

In  this  study,  the  location  of  the  gas-liquid  interface  along  the 
flow  channel  direction  at  various  cathode  humidity  conditions 
and  its  effect  on  cell  performance  were  elucidated  by  modeling  a 
3D  PEMFC  system  using  CFDRC.  Fig.  2  shows  the  effect  of  the 
relative  humidity  of  the  cathode  on  the  location  of  the  interface 
where  the  liquid  water  begins  to  condense  along  the  flow  channel 
at  a  cell  operating  voltage  of  0.7  V.  The  interface  is  defined  as  the 
location  where  liquid  water  begins  to  condense.  The  horizontal 
dotted  line  indicates  the  interface  between  the  flow  channel  and 
the  gas  diffusion  layer.  A  higher  cathode  relative  humidity  cor¬ 
responds  to  a  smaller  distance  between  the  gas-liquid  interface 
and  the  gas  inlet,  as  clearly  displayed  in  Fig.  2(a).  Furthermore, 
the  gas-liquid  interface  moves  to  the  flow  channel  inlet  as  the 
relative  humidity  of  the  cathode  increases.  This  phenomenon  is 
caused  by:  (i)  the  decrease  in  the  amount  of  evaporated  water 
through  the  flowing  gas  stream;  (ii)  the  increase  in  the  partial 
pressure  of  water  and  the  ability  to  reach  the  saturation  pressure 
of  the  vapor  water  relatively  quickly  to  form  liquid  water  earlier 
as  the  relative  humidity  of  the  cathode  increases.  On  the  contrary, 
the  gas-liquid  interface  moves  closer  to  the  catalyst  layer  when 
the  relative  humidity  of  the  cathode  is  less  than  60%,  as  shown  in 
Fig.  2(b).  In  conclusion,  increasing  the  relative  humidity  of  the 
cathode  can  reposition  the  gas-liquid  interface  and  cause  liquid 
water  to  appear,  affecting  the  performance  of  the  cell,  because 


(b)  Z  ( m  ) 


dw  dCk 

u  =  —  =  v  =  -  =  0 

3  z  d  z 


(24) 


Fig.  2.  Effects  of  relative  humidity  of  cathode  on  the  location  of  the  interface 
where  liquid  water  begins  to  condense  along  the  flow  channel  at  a  cell  operating 
voltage  of  0.7  V.  (a)  RHca  =  20-100%  and  (b)  RHca  =  20-60%. 
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Fig.  3.  Effect  of  relative  humidity  of  the  cathode  on  cell  performance  at  cell 
operating  voltage  of  0.7  V.  (a)  I-V  curves  and  (b)  I-P  curves. 

liquid  water  may  occupy  the  pores  in  the  porous  media,  reducing 
the  amount  of  fuel  gas  that  can  reach  the  catalyst  layer. 

Fig.  3(a)  and  (b)  plot  I-V  and  I-P  curves  at  various  relative 
humidities  of  the  cathode  at  an  operating  voltage  of  0.7  V.  The 
results  in  Fig.  3(a)  reveal  that  the  cell  performance  decreases 
as  the  relative  humidity  of  the  cathode  increases,  because  the 
amount  of  liquid  water  increases  with  the  relative  humidity. 
Accordingly,  the  pores  in  the  porous  media  were  obstructed  by 
liquid  water  at  the  cathode-side,  reducing  the  amount  of  reaction 
gas  to  the  catalyst  layer.  Therefore,  the  performance  of  the  cell 
gradually  decreases.  Also,  the  polarization  curves  do  not  seem 
to  vary  as  the  relative  humidity  of  the  cathode  increases,  at  an 
operating  voltage  that  exceeds  ~0.65  V.  A  large  voltage  results 
in  a  small  current  density,  and  therefore,  a  relatively  small  elec¬ 
trochemical  reaction  rate.  However,  the  performance  of  the  cell 
decreases  significantly  as  the  operating  voltage  declines  below 
the  ^0.65  V  operating  voltage,  since  a  large  current  density 
accelerates  the  rate  of  the  electrochemical  reaction.  In  particular, 
in  Fig.  3(a),  the  higher  cathode  relative  humidity  (100%)  offers  a 
better  cell  performance  than  the  lower  one  (20%)  at  an  operating 
voltage  of  over  ~0.65  V,  because  the  rate  of  the  electrochemical 
reaction  decreases  as  the  operating  voltage  increases.  There¬ 
fore,  the  water  content  of  the  membrane  in  the  initial  stage 
increases  with  the  relative  humidity.  Thus,  the  cell  performance 
is  improved  as  the  relative  humidity  of  the  cathode  increases. 


(a)  Z  (  m  ) 


Fig.  4.  Effect  of  cell  operating  voltage  on  the  location  of  the  interface  where 
liquid  water  begins  to  condense  along  the  flow  channel  at  a  relative  humidity  of 
the  cathode  of  80%.  (a)  V=  0.5-0. 9;  (b)  V  =  0. 8-0.9. 

By  contrast,  at  operating  voltages  of  below  ~0.65  V,  the  cell 
performance  improves  as  the  relative  humidity  of  the  cathode 
decreases,  because  accelerating  the  electrochemical  reaction 
produces  more  water,  increasing  the  water  content  in  the  mem¬ 
brane  at  a  low  cathode  relative  humidity,  such  as  20%.  However, 
at  a  high  relative  humidity,  flooding  may  occur  on  the  cathode- 
side.  Additionally,  the  power  density  increases  as  the  relative 
humidity  of  the  cathode  decreases,  as  shown  in  Fig.  3(b). 

Fig.  4  shows  the  effects  of  various  operating  voltages  on  the 
location  of  the  gas-liquid  interface  along  the  flow  channel  at  a 
relative  humidity  of  the  cathode  of  80%.  Fig.  4  reveals  that  the 
gas-liquid  interface  is  close  to  the  catalyst  layer  at  a  high  oper¬ 
ating  voltage,  indicating  that  less  water  is  generated  at  a  lower 
electrochemical  reaction  rate.  However,  since  more  water  is  gen¬ 
erated  at  a  high  electrochemical  reaction  rate,  the  gas-liquid 
interface  is  close  to  the  gas  diffusion  layer  and  the  flow  channel  at 
a  low  operating  voltage.  In  closing,  the  gas-liquid  interface  grad¬ 
ually  moves  to  the  gas  inlet  as  the  operating  voltage  decreases, 
because  reducing  the  operating  voltage  increases  the  current 
density.  Therefore,  the  electrochemical  reaction  rate  increases. 

Figs.  5  and  6  plot  the  oxygen  and  water  fractions  in  the  gas 
flow  channel  and  the  gas  diffusion  layer,  respectively,  of  the 
cathode-side  in  the  direction  of  the  flow  channel  at  an  oper¬ 
ating  voltage  of  0.7  V  and  relative  humidity  of  the  cathode  at 
80%.  Larger  oxygen  fraction  appears  in  the  gas  inlet,  and  then 
gradually  decreases  in  the  flow  direction,  as  shown  in  Fig.  5.  The 
results  reveal  that  the  decrease  in  the  oxygen  fraction  is  caused  by 
the  catalyst  layer  consuming  more  oxygen  on  the  cathode-side. 
Hence,  the  oxygen  fraction  is  lowest  near  the  catalyst  layer.  By 
contrast,  Fig.  6  reveals  that  the  water  fraction  gradually  increases 
along  the  flow  direction,  because  of  the  production  of  water  in  the 


Chun-I  Lee,  H.-S.  Chu  /  Journal  of  Power  Sources  161  (2006)  949-956 


955 


Fig.  5.  Oxygen  fraction  in  the  cathode  gas  channel  and  gas  diffusion  layer  along 
the  flow  channel  at  a  cell  voltage  of  0.7  V  and  a  relative  humidity  of  the  cathode 
of  80%. 


Fig.  6.  Water  fraction  in  the  cathode  gas  channel  and  the  gas  diffusion  layer 
along  the  flow  channel  at  a  cell  voltage  of  0.7  V  and  a  relative  humidity  of  the 
cathode  of  80%. 

catalyst  layer  of  the  cathode  by  the  electrochemical  reaction,  as 
well  as  the  transport  of  water  from  the  anode- side  to  the  cathode- 
side  by  electro-osmotic  drag.  Accordingly,  the  water  fraction  is 
the  highest  near  the  catalyst  layer.  Since  the  water  mass  fraction 
is  higher  near  the  catalyst  layer,  where  the  corresponding  water 
partial  pressure  exceeds  the  saturated  vapor  pressure,  causing 
the  formation  of  liquid  water.  Therefore,  Fig.  7  shows  the  liquid 
water  saturation  field  in  the  gas  channel  and  the  gas  diffusion 
layer  of  the  cathode- side  at  an  operating  voltage  of  0.7  V  and  a 
relative  humidity  of  the  cathode  of  80%.  The  figure  shows  that 
saturation  of  the  liquid  water  increases  along  the  flow  channel, 
because  the  electrochemical  reaction  causes  the  partial  pressure 
of  the  water  to  exceed  the  saturated  vapor  pressure,  causing  liq¬ 
uid  water  to  condense.  The  capillary  force  also  causes  liquid 


Fig.  7.  Liquid  water  saturation  field  in  the  cathode  gas  channel  and  gas  diffusion 
layer  at  a  cell  voltage  of  0.7  V  and  a  relative  humidity  of  the  cathode  of  80%. 


water  to  move  toward  the  gas  diffusion  layer.  Hence,  this  inter¬ 
face  is  the  single-phase  region  when  the  saturation  of  the  liquid 
water  is  0.  It  becomes  a  two-phase  region  when  the  saturation  of 
the  liquid  water  exceeds  0.  According  to  Fig.  2,  to  the  right  of 
the  point  where  water  condensation  starts,  a  two-phase  region  is 
present  along  the  flow  channel.  A  single-phase  region  exists  on 
the  other  side. 

4.  Conclusions 

A  multi-dimensional,  multi-component,  computational  fluid 
dynamic  model  was  developed  to  study  the  effect  of  the  cath¬ 
ode  humidification  conditions  on  the  location  of  the  interface 
where  the  liquid  water  begins  to  condense  along  the  flow  chan¬ 
nel  in  a  PEM  fuel  cell.  The  model  results  support  the  following 
conclusions. 

(1)  The  gas-liquid  interface  moves  from  the  catalyst  layer 
toward  the  diffusion  layer  as  the  relative  humidity  of  the 
cathode  increases.  When  the  relative  humidity  of  the  cath¬ 
ode  reaches  100%,  the  gas-liquid  interface  is  close  to  the 
gas  flow  channel  inlet. 

(2)  As  the  condensed  water  in  the  pores  in  the  porous  media 
blocks  the  transport  of  fuel  gas,  the  cell  performance  and 
power  density  decrease  as  the  relative  humidity  of  the  cath¬ 
ode  increases. 

(3)  Reducing  the  cell  operating  voltage  reduces  the  distance 
between  the  gas-liquid  interface  and  the  gas  flow  channel 
inlet,  because  the  higher  current  is  associated  with  a  higher 
electrochemical  reaction  rate. 

(4)  The  decreased  oxygen  fraction  and  the  increased  water  frac¬ 
tion  along  the  flow  channel  are  related  to  the  electrochemical 
reaction  of  the  cathode  catalyst  layer,  as  oxygen  is  consumed 
and  water  produced. 

(5)  The  transport  of  liquid  water  through  the  porous  media  is 
driven  by  the  shear  force  of  gas  flow  and  the  capillary  force. 
Accordingly,  liquid  water  moves  from  the  catalyst  layer 
toward  the  gas  diffusion  layer. 
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